In-plane magnetodrag in dilute bilayer two-dimensional systems: a Fermi liquid theory 
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Motivated by recent experimental results reporting anomalous drag resistance behavior in dilute 
' bilayer two-dimensional (2D) hole systems in the presence of a magnetic field parallel to the 2D plane, 

' we have carried out a many-body Fermi liquid theory calculation of bilayer magnetodrag comparing 

, it to the corresponding single layer magnetoresistance. In qualitative agreement with experiment we 

04 ' find relatively similar behavior in our calculated magnetodrag and magnetoresistance arising from 

^jr^, the physical effects of screening being similarly modified ("suppressed") by carrier spin polarization 

1-^ . (at "low" field) and the conductivity effective mass being similarly modified ("enhanced") by strong 

' magneto-orbital correction (at "high" fields) in both cases. We critically discuss agreement and 

, disagreement between our theory and the experimental results, concluding that the magnetodrag 

' data are qualitatively consistent with the Fermi liquid theory. 

PACS Number : 73.40.-c, 73.21.Ac, 73.40.Kp 
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Much attention has recently focused on low-density 2D 
systems in semiconductor structures [1-7] where carrier 
transport properties may be strongly affected by inter- 
action effects. In particular, low temperature transport 
[1,2], magnetotransport [3,4], drag [5,6], and magneto- 
drag [7] properties have recently been studied in low- 
density electron [1,3,6] and hole [2,4,5,7], single-layer 
[1-4] and bilayer [5-7] systems, providing a great deal 
of detailed quantitative information on the tempera- 
ture, density, and magnetic field dependent 2D resistivity 
p{T,n,B) and 2D drag-resistivity pD{T,n, B) behavior. 
(The externally applied magnetic field being discussed 
throughout this work and in the relevant experimental 
references [1-7] is an 'in-plane' magnetic field B applied 
parallel to the 2D layer.) A very recent experimental 
work [7] by Pillarisetty et. al. reports some striking qual- 
itative resemblance between the 2D bilayer drag pu and 
the corresponding single-layer resistivity p as a function 
of the applied parallel field S in a low density low disor- 
der 2D GaAs hole system. Since the physical mechanisms 
underlying p and pu are generally thought to be quali- 
tatively different at low temperatures, the experimental 
observations of ref. 7 take on important qualitative sig- 
nificance. In particular, p in high-mobility 2D systems at 
low temperatures is entirely due to scattering by random 
charged impurity centers whereas pu at low temperatures 
arises entirely from inter-layer electron-electron scatter- 
ing. (Electron-phonon scattering makes negligible contri- 
butions to both p and po at low temperature [1,2].) Since 
electron-electron scattering does not directly contribute 
to p in translationally invariant 2D semiconductor sys- 
tems, the reported [7] qualitative similarity between the 
observed p and po behaviors in its magnetic field depen- 
dence presents a significant theoretical challenge. Since 
it is manifestly obvious that electron-impurity scatter- 
ing can at best play an unimportant and indirect sec- 
ondary role [8] in determining the interlayer drag resis- 
tance, the experimental observation of ref. 7 raises very 
serious fundamental questions regarding our understand- 



ing of the nature of the ground state of a low-density 2D 
carrier system. We note that electron-electron interac- 
tion induced umklapp scattering, which could contribute 
to the single-layer resistivity (since umklapp processes do 
not conserve momentum), is completely irrelevant in 2D 
semiconductor structures where all the electronic physics 
occurs essentially at the zone-center F point in the ef- 
fective mass approximation sense (and the real lattice 
structure does not play any role) . 

In view of the considerable fundamental significance 
of the issues raised by the experimental observations, we 
present in this Letter a careful theoretical calculation of 
both p{B) and pd{B) in a low-density 2D carrier system 
within the canonical many-body Fermi liquid theory that 
has earlier been found to be successful in providing a rea- 
sonable qualitative (and perhaps even semi-quantitative) 
description of the approximate temperature and density 
dependence of p [9] and po [10] at low temperatures and 
densities in the absence of any applied in-plane magnetic 
field. We note that the zero-field temperature and carrier 
density dependence of 2D resistivity p(T, n) and 2D drag 
resistivity pD{T,n) in the absence of any external mag- 
netic field are certainly very different as one would ex- 
pect on the basis of p and pD being determined by differ- 
ent scattering processes: p by screened charged impurity 
scattering and pu by interlayer electron-electron scatter- 
ing. For example, p shows [1-4] an approximate linear 
increase with T at low temperatures as is expected [9] 
for screened Coulomb impurity scattering and po shows 
[5,8,10] an approximate quadratic increase with T at low 
temperatures as is expected for electron-electron scatter- 
ing. (Similarly the carrier density dependence of p and 
PD are also very different at S = 0.) The question there- 
fore naturally arises why the in-plane magnetic field de- 
pendences of p{B) and pd{B) reported in ref. [7] show 
qualitative similarities. 

We theoretically argue, showing concrete calculated re- 
sult within the many body Fermi liquid theory, that the 
qualitative magnetic field dependence of p{B) and pd{B) 
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should indeed be similar since the scattering processes 
controlling the two properties (electron charged impu- 
rity scattering for p and electron-electron scattering for 
Pd) arc both screened by the carriers themselves and the 
dominant behavior in both cases arises primarily from 
the magnetic field dependence of electronic screening [11] 
(through the spin polarization process) and (somewhat 
to a lesser degree) from the magneto-orbital effect [12] 
(through the modifications of the 2D conductivity ef- 
fective mass and the confining quasi-2D wave function). 
The reported qualitative similarity between magnetodrag 
and magnetoresistance thus arises from drag and resis- 
tance being dominated by screened carrier-carrier scat- 
tering and screened carrier-impurity scattering respec- 
tively. The fact that long-ranged charged impurity poten- 
tial is the dominant source of resistive scattering in 2D 
semiconductor structures (and this long-ranged charged 
impurity scattering must necessarily be screened by the 
carriers) is therefore the key reason for the broad qual- 
itative similarity between pd{B) and p{B) reported in 
[7]. 

We start by writing down the zero-field theoretical for- 
mulae for p [9] and po [10,8] in the many-body Fermi liq- 
uid RPA-Boltzmann theory approximation widely used 
in the literature. The resistivity is given by p~^ = 
ne'^{T)/m, where n, m are the 2D carrier density and 
the conductivity effective mass respectively whereas the 
transport relaxation time r is given by 

= "^^ritlueiOi - k')f{l - coseu.i,')S{ek - Sk'), 

(1) 

with (r) being a thermal average over the carrier energy 
e. Here rij is the density of charged impurity centers in 
the 2D system (including the interface and the insula- 
tor), and Uei{q) is the screened carrier- impurity scatter- 
ing strength give by u{q) — v'^{q)/e{q), where e{q) is the 
single-layer 2D carrier dielectric function. (For details on 
the derivation and implications of the formula for p, see 
ref. [9].) The drag resistivity is given by 

fq^d^q r dm F,iq,cu)F2{q,L0) 
27re^n^kBT J {27rf J 27t sinh2(/3a;/2) ' ^ ' 

where ^"1^2(9,^) = |mi2(<?, ^)|Iinnii_22(q, t^), with uf^ = 
■^1122/^(91'^) is the dynamically screened interlayer 
Coulomb interaction between layers 1 and 2, and 11 is 
the 2D polarizability. (We consider the so-called bal- 
anced situation here with the same carrier density n in 
both layers.) Note that the dielectric function €{q, w) = 
1 — v{q)Il{q,Lu) entering Eq. (2) is the two component 
dielectric tensor for the bilayer system [13]. (For details 
on the drag formula and its implications, see refs. [8,10].) 

It is important to emphasize that dielectric screening 
by the carriers themselves is a key ingredient in deter- 
mining both p and po although the static single layer 



(scalar) dielectric function e(q) determines p through the 
screened charged impurity scattering whereas the dynam- 
ical bilayer (tensor) dielectric function e{q, u) determines 
Pd through the screened interlayer Coulomb interaction. 
At low carrier densities used in ref. [7], the difference 
between static and dynamic screening is not of any qual- 
itative significance since the effective plasma frequency 
scale is rather low at low densities. Therefore both p and 
pn depend on the carrier dielectric function properties, 
which is why they have qualitative similar magnetic field 
dependence as we show and discuss below. 

We have carried out a thoroughly nontrivial general- 
ization of the above theories for p and po to the finite in- 
plane magnetic field situation p{B), pd{B). Details will 
be provided elsewhere [14], but here we mention the main 
physical effects of the applied field for p and po- The ap- 
plied field has two completely different physical effects 
through its coupling to carrier spin ("magneto-spin") 
[11] and orbital ("magneto-orbital") dynamics [12]. The 
magneto-spin effect arises from field-induced carrier spin 
polarization due to the Zeeman coupling, and saturates 
at a density dependent saturation field Bs when the car- 
rier system is fully spin polarized (i.e. the magneto-spin 
effect exists only for B < Bs). The magneto-orbital ef- 
fect [12] arises from the orbital coupling of the in-plane 
magnetic field to the transverse dimension due to the 
quasi-2D nature of the 2D layer and the magneto-orbital 
effect is therefore monotonically increasing with increas- 
ing magnetic field since this orbital coupling is important 
only when the magnetic length / = {ch/ eB)^^"^ is smaller 
than the quasi-2D width of the 2D system. Thus one im- 
portant qualitative difference between the magneto-spin 
and the orbital effect is that the spin effect is essentially a 
"weak-field" effect lasting only upto the saturation field 
Bs whereas the magneto-orbital effect increases mono- 
tonically with increasing field. 

The magneto-spin mechanism itself has two distinct ef- 
fects: Suppression of screening due to spin polarization 
[11] as the spin degeneracy decreases from 2 (at _B = 0) to 
1 (at £? > Bs) and the increase of the effective 2D Fermi 
surface as the value of the 2D Fermi wave vector kp in- 
creases by a factor of with B increasing from zero to 
Bs due to the lifting of the spin degeneracy. Similarly, 
the magneto-orbital mechanism also has two distinct ef- 
fects: The increase of the transport effective mass in the 
direction perpendicular to the magnetic field direction 
and the field-induced intersubband scattering among the 
quasi-2D subband — both of these are only operational 
at relatively high fields when I < a where a is the aver- 
age transverse width of the carrier wave function. It is 
important to realize that three of these four field induced 
effects (spin polarization induced screening suppression, 
and both of the magneto-orbital effects) always produce 
positive magnetoresistance whereas the Fermi surface ef- 
fect (which is significant only at high carrier densities 
where 2kF ^ qrp, qTF being the screening wave vector) 
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1. Calculated (a) magnetoresistance p[B\\) and (b) 



magnetodrag pd(S||) for hole density p = 2.15 x 10 cm at 
various temperatures as a function of parallel field B\\ . 



always produces a negative magnetoresistance. For the 
hole-doped low-density samples of rcf . [7] , the Fermi sur- 
face (i.e. kp V^kp as B ^ Bs) effect is negligible 
since the system is in the strong screening qtf ^ ^kp 
limit. 

The combination and the interplay of magneto-spin 
and magneto-orbital effects are quite complex and sensi- 
tive to the parameter {n,T. B) details, but a few general 
comments can still be made: (1) At low carrier densities 
[7] of interest to us, the static and dynamic screening 
operational respectively in p{B) and pd{B) behave simi- 
larly, and therefore the spin-polarization induced screen- 
ing effect is qualitatively similar for p{B) and Pd{B); (2) 
since field-induced magneto-spin effect operates only for 
B < Bg, both p{B) and pd{B) manifest a cusp- type 
structure at B = Bg where spins are completely po- 
larized; (3) the maximum theoretically allowed magne- 
toresistance p{B)/p{0) and magnetodrag pd{B)/pd{o) 
corrections arising from the spin polarization induced 
magneto-screening mechanism arc factors of 4 and 16 re- 
spectively since screening itself could be suppressed at 
most by a factor of 2 due to spin polarization effect (and 
Eqs. 1 and 2 respectively for p and po come with the 
second and the fourth power of the spin degeneracy) ; (4) 
the main magneto-orbital effect for the relatively narrow 
p-GaAs quantum well systems (width ^ 150 A) used in 
ref. [7] is the enhanced conductivity mass at higher mag- 
netic field values — the condition I ^ 150 A necessary 
for strong magneto-orbital correction is satisfied for B ::g> 
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FIG. 2. Calculated (a) magnetoresistance p(S||) and (b) 
magnetodrag pd{B\\) at a temperature T = 100 mK for dif- 
ferent density p =1.2, 1.5,1.75, 2. 15x 10^° cm~^ as a function 
of parallel field Bu . 



4 T whereas the spin polarization saturation field Bs for 
the low hole densities used in ref. [7] is ~ 3 — 6 
T; thus the magneto-spin effects dominate for B upto 
3 — 6 T whereas the magneto-orbital effects dominate at 
higher fields; (5) the magneto-orbital effects are "similar" 
in both cases since both p and pn are proportional to the 
field-dependent effective mass (which increases quadrat- 
ically with the applied field). We mention that these 
five features are in excellent qualitative agreement with 
experimental results [7] . 

In Figs. 1 - 4 we show our calculated results for p 
and Pn within the RPA-Boltzmann Fermi liquid theory. 
Our theory incorporates all realistic effects [9-11] with 
the charged impurity density (n^) determining p as the 
only unknown free parameter. Our results in Figs. 1 and 
2, where p{B) and po{B) are shown for different temper- 
atures and different densities, respectively, bear excellent 
qualitative resemblance to the corresponding experimen- 
tal results in ref. [7]. We are not claiming quantitative 
agreement with experiment by any means since our the- 
ory is necessarily approximate at the low carrier densi- 
ties used in ref. [7] since no exact description of corre- 
lation effects at low densities exists for interacting quan- 
tum Coulomb system of interest here. The qualitative 
agreement between theory and experiment is, however, 
obviously apparent even on a casual comparison between 
our Figs. 1 and 2 and the corresponding Figs. 1 and 
2 in ref. [7]. In particular, both theory and experiment 
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FIG. 3. Calculated (a) p{B) and (b) Pd{B) including only 
the magneto-spin or only the magneto-orbital effects for a hole 
density p = 2.15 x 10^°cm-^ at T = WOmK. 

manifest qualitatively similar, but by no means identical, 

behaviors in p{B) and po{B), arising, as argued above, 
from the magneto-spin and magneto-orbital effects. 

It is worthwhile to theoretically consider the orbital 
and the spin effects separately. (Experimentally this can- 
not, of course, be done but one could get some approx- 
imate idea about the relative behavior of the magneto- 
spin and the magneto-orbital effects in p{B) and pij{B) 
by concentrating on the 'low' B (< Bg) and the 'high' 
B (> Bg) regimes, respectively.) In Fig. 3 we show the 
calculated p{B) and pd{B) including only the magneto- 
spin or only the magneto-orbital effects. Again the im- 
portance of the 'low-field' magneto-spin and the 'high- 
field' magneto-orbital effects on both p{B) and Pd{B) 
are manifestly obvious in our theoretical results. 

Finally in Fig. 4 we present some clear-cut theoreti- 
cal predictions for the temperature dependence of bilayer 
magnetodrag pd{B; n, T) in the presence of the in-plane 
magnetic field B. In particular, we fit the temperature- 
dependence of pd at a fixed low density to approximate 
power law behaviors: Pd{T) ~ with the magnetic 
field dependent exponents a{B) indicating the nature of 
the temperature dependence of magnetodrag at various 
magnetic field values. The striking theoretical prediction, 
which stands out in Fig. 4, is that a manifests a very 
strong magnetic field dependence with a{B) decreasing 
from a low-field value of about 2.3 (for B < Bg) to a 
high-field value of about 1.8 (for B > Bg). This sharp 
drop in the temperature exponent of magneto-drag is a 
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FIG. 4. a vs. S|| for hole density p = 2.15 x 10^° cm"^, 
where the exponent a{B^^) is deduced from the linear fit of 
calculated magnetodrag, pD(B||,r) = cT". 



direct consequence of the strong suppression in magneto- 
screening arising from the carrier spin polarization in- 
duced by the in-plane magnetic field. (Note that these 
exponents are 'effective' exponents and not exact expo- 
nents.) Our approximate analysis of the experimental 
data [7] indicate that our theoretical results for a shown 
in Fig. 4 are in excellent qualitative (and reasonable 
quantitative) agreement with rcf. [7], where a. changes 
from around 2.5 for S ~ to about 1.3 for large B 

We conclude by emphasizing that our Fermi liquid the- 
ory based detailed calculations are in excellent qualitative 
agreement with the experimentally observed magneto- 
drag data [7], and therefore more exotic non- Fermi liq- 
uid theory [15] descriptions (which cannot typically pro- 
duce quantitative results as shown in our Figs. 1 — 4) 
seem unnecessary. The 'smoking gun' breakdown of the 
Fermi liquid description of bilayer drag experiment would 
be the observation of a drag resistance which remains 
finite as T — > since within the Fermi liquid theory 
Pd(T = 0) = 0. All existing experimental data seem 
to be consistent with the Fermi liquid conclusion that 
poiT ^ 0) ^ 0. 

This work is supported by the US-ONR, the NSF-ECS, 
the ARO, the ARDA, and the LPS. 
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